Column collapse of granular rods 
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We investigate the collapse of granular rodpiles as a function of particle (length/diameter) and 
pile (height/radius) aspect ratio. We find that, for all particle aspect ratios below 24, there exists a 
critical height Hi below which the pile never collapses, maintaining its initial shape as a solid, and 
a second height Hu above which the pile always collapses. Intermediate heights between Hi and 
Hu collapse with a probability that increases linearly with increasing height. The linear increase in 
probability is independent of particle length, width, or aspect ratio. When piles collapse, the runoff 
scales as a piecewise power-law with pile height, with rf ~ jc/i-2±o.i ^^^ p-jg }ieights below He ~ 0.74 
and rf « J/O-Sioi for taller piles. 



PACS numbers: 45.70.Ht 



I. INTRODUCTION 

The collapse and flow of granular materials is funda- 
mental to catastrophic events such as avalanches, land- 
slides, and pyroclastic flows. Understanding how initially 
stationary piles of granular material collapse, how far the 
material travels before stopping, and the final pile ge- 
ometry, is an important step in efforts to deal with the 
practical consequences of such natural phenomena. 

Previous experimental investigations in granular col- 
umn collapse [l|, 01 have focused on cylindrical piles (ini- 
tial radius ro) of varying height H . The radius of the col- 
lapsed pile, r/, scaled as f = (r/ ~r'o)/ro, shows a power- 
law dependence on dimensionless height H = H/ro- 
There is a crossover in the power-law exponent when H 
becomes larger than a critical height He, although there 
is not consensus on the specific value of He- Lube et al. 
P,[l found 



H^- H<1.7 
H^-^ H > 1.7 



(1) 



for a variety of irregular granular materials, independent 
of surface roughness. The crossover was also observed to 
coincide with a slight change in the final pile shape. 

Lajeunesse et al. [3, 3 studied more regular mate- 
rials (smooth glass beads) and also observed a change 
in scaling of the runoff, but found the critical height to 
be He ~ 0.74. The collapse of low {H < 0.74) piles 
resulted in a mesa-like final state, with some fraction 
of the top surface remaining horizontal. The horizon- 
tal plateau shrank with increasing H, finally disappear- 
ing for H > 0.74. Both Lube and Lajeunesse found 
little change in behavior for a variety of particle types 
or surface roughness. Their studies encompassed fine 
and coarse quartz sand, salt, couscous, rice, sugar, and 
smooth glass beads. 



It is well known that rod-like granular materials of suf- 
ficient aspect ratio a = L/d can support 90° angles of re- 
pose, and thus piles of these granular materials maintain 
the shape of the container in which they were initially 
formed (c.f. figure [D (right)). Philipse [^-Q first studied 
this solid plug-like behavior, and observed that particles 
with aspect ratio greater than « 35 could not be poured 
from a container, but rather emerged as a solid, a behav- 
ior has been assumed to be due to geometric entangle- 
ment of the particles. 




1svfspsOrit.edu| |http : //piggy • rlt ■ edu/f ranklln/ 1 



FIG. 1: Rod-like piles of moderate height may collapse (left), 
but may also retain their original, cylindrical shape (right). 
Shown are results from two experiments on aspect ratio 
a = 16 particles with identical initial heights. We find the 
probability of collapse increases linearly with pile height. 

This has led to a great deal of research on rods of 
various aspect ratio. Villarruel et al. Q found that 
rods of aspect ratio sa 4 constrained to a cylinder would, 
when tapped, align with the cylinder walls and compact. 
Lumay and Vandewalle [l3| extended this experiment to 
much larger aspect ratios, and developed a 2d lattice 
model that reproduced the salient features of the com- 
paction curves. Lumay and Vandewalle subsequently [ll| 
used a two-dimensional experiment to highlight the im- 
portance of rotation in rod rearrangements. Blouwolff 
and Fraden [l^] investigated the contact number in 3d 
packings, and Stokely et al. [iJl looked at static proper- 
ties of 2d packings. Pournin et al. [IJ, [l^ have written 
discrete element method (DEM) simulations of sphero- 
cylindrical particles that show crystallization and order- 
ing under various excitations. Recently, Hidalgo et al. 
[16] looked at the role particle shape plays on stress prop- 
agation through a granular packing, while Azema and 
Radji[l3| looked at stress-strain relations of rod pack- 
ings under shear. Finally, there has been some attention 



paid to the nature of random packings. Wouterse et al. 
|18| looked at the microstructure of random packings of 
spherocyhnders and Zeravcic et al. [l^ have looked at 
the excitations of random ellipsoid packings at the onsite 
of jamming. Despite all this attention, there have been 
no structural characteristics identified that explain the 
increased rigidity of large aspect-ratio particles. 

Desmond et al. [20[ studied the transition from clas- 
sically granular to solid-like behavior by measuring the 
force on a small object pulled through a rodpile. They 
found characteristically granular slip-stick behavior for 
low aspect ratio particles, while large aspect ratio par- 
ticles responded as a solid, moving as a single coher- 
ent unit. Desmond et al. were the first to observe the 
existence of a transition region separating the granular 
and solid behaviors. Moderate aspect-ratio particles dis- 
played both granular and solid behavior, often alternat- 
ing between the two. Typically a pile would briefly dis- 
play solid-like behavior and move coherently, before col- 
lapsing around the intruder. Desmond et al. mapped the 
behavior in a two-dimensional parameter space of parti- 
cle aspect ratio and inverse container size, but there was 
no quantitative measure of the transition region. 



ro)/rQ. As shown in Fig. [2l sand, acrylic rods of aspect 
ratio 8, 10, and 16, and Teflon rods of aspect ratio 8 and 
16 all show a power-law dependence on normalized pile 
height H = H/tq with change in exponent at a critical 
height He- To determine the the exponents and crossover 
height we perform a 3-parameter least-squares fit of the 
data to the piece- wise power-law function 



i7"i H < H^ 



(2) 



The minimization is done analytically for ai and a2 and 
numerically for He- We have done this fitting both for 
every individual data run, averaging the resulting expo- 
nents and transition point, as well as for the aggregate 
collection of all data. Both analyses fall within the same 
standard-deviation uncertainty for the parameters, with 
ai = 1.2 ± 0.1, aa = 0.6 ± 0.1, and H, = 1-1 ± 0.3, or, 
putting these values into Eq. [2] 



jji.2±o.i jj ^ 1.1 ±0.3 
^o.6±o.i ^ > 1.1 ±0.3 



(3) 
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II. EXPERIMENTAL SETUP 

Our experimental apparatus follows that of Lube and 
Lajeunesse [l|-lj]. An acrylic cylinder with inner ra- 
dius To is filled with material to height H and then 
quickly lifted vertically. Two cylinder radii were tested, 
rg — 5.4 ± 0.1 cm and 7.3 ± 0.1 cm. We observed no de- 
pendence on cylinder radius and, unless otherwise noted, 
all reported results are from the smaller cylinder. The 
cylinder was filled by dropping particles in from the top. 
Two orthogonal cross-beams spanned the top entrance to 
the cylinder and randomized the particles' orientation as 
they fell into the cylinder. The cross-beams also broke 
up any clumps that existed prior to pouring. 

The cylinder was raised with speeds ranging from 1-10 
cm/s; little dependence on cylinder speed was observed. 
After the pile had collapsed, the distance from the cylin- 
der's central axis to the pile edge, defined as the midpoint 
of the farthest particle in contact with the main pile, was 
measured in four directions (ix, ±y) and averaged. 

Rods were cut from acrylic and Tcfion rod stock. Par- 
ticle diameters were D — 0.16 ± 0.01 — 0.95± cm, lengths 
L — 2. 5 ±0.2 — 7. 6 ±0.2 cm, resulting in a range of aspect 
ratio from 2.6 — 47.5. The nominal coefficient of friction 
acrylic is fiacr = 0.8; that of Tefion jit = 0.04. Ordi- 
nary sand was tested as a control, and compared with 
the results of Lube and Lajeunesse et al. 



III. SCALING OF RUNOFF WITH HEIGHT 

Following Lube, we define the runoff as the change in 
pile radius, normalized by the initial radius: f = (r/ — 
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Normalized pile height (H = H/r ) 



FIG. 2: The runoff distance f of a collapsed granular pile 
shows a power-law dependence on initial height H, with a 
transition between linear and square-root scaling. Both ra- 
dius and pile height H are made dimensionless by the initial 
pile radius tq. For H <« 1.1 the runoff increases linearly 
with H, whereas for larger piles the runoff increases as H^'^ 
(solid lines are power law fits). Both scaling behaviors, as well 
as the transition point, arise from conservation of volume ar- 
guments that assume all material outside an internal, stable 
cone collapses. 



The range of heights tested is only about a decade for 
the lower region, and less than a decade for the upper 
region, and so we only assert that our results are consis- 
tent with the findings of Lube et al. [2|| and Lajeunesse 
et al 0,13. 

As worked out by Lajeunesse Q, all of the behavior in 
Fig. [5] — the two power-law scaling exponents and the 



crossover between the two — can be explained by conser- 
vation of volume, assuming the material does not com- 
pact as it collapses. The crossover height occurs when 
the dimensionless pile height is such that 

Hq = tan 6'c, 

where 9c is the dynamic angle of repose. Our observed 
transition height of Hq — 1.1 ± 0.3 implies an angle of 
repose of about 9c ~ 47° ± 7°, slightly higher than that 
found by Lajeunesse for round and irregular particles, 
indicating a slightly increased structural stability of rod- 
piles. 

To conclude this section, the runoff of rods of even 
moderate aspect ratios follows similar scaling as that of 
sand. The comparable transition point indicates a sim- 
ilar angle of repose suggesting that, when the pile col- 
lapses, the cylindrical particles roll much like spheres. 
The square-root scaling of runoff at large pile heights in- 
dicates that the dynamics of pile collapse are also similar, 
with an immobile inner core setting the final pile height. 



IV. SOLIDITY OF MODERATE 
ASPECT-RATIO PARTICLES 



exist to guide the eye and reveal qualitative changes in 
behavior when the pile height is less than L/4, between 
L/4 and 3L/4, and greater than 3L/4. It should be noted 
that the data in Fig. [3] come from a variety of particle 
diameters and initial cylinder radii, and so argue that it 
is the particle length that determines the maximum pile 
height. Reducing particle friction narrowed, but did not 
eliminate, the transition region. 
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One of the more striking behaviors of large aspect-ratio 
granular material is their ability to maintain the shape 
of their initial container. As expected, piles formed from 
particles with aspect ratio (L/d) greater than 24 never 
collapsed, maintaining their cylindrical shape. While 24 
is slightly lower than the value commonly assumed nec- 
essary to maintain plug- like behavior (0, 0, HI, [201), it 
is consistent with the finding of Desmond et al. that 
finite-size container effects can induce plug-like behav- 
ior at lower aspect ratios [20|. Surprisingly, however, we 
find that this solid-like behavior manifests in moderate 
aspect-ratio particles as well (see Fig. [T] (right)). 

The behavior of piles formed from particles with inter- 
mediate aspect ratios 4 < a < 24 could not be predicted 
a priori, with experiments with visually equivalent ini- 
tial conditions producing dramatically different results. 
Shown in Fig. [T] are two piles that resulted with parti- 
cles of aspect ratio a = 16 with relatively low initial pile 
heights. As can be seen, one pile collapsed, resulting in 
a mound formation with a relatively well-defined angle 
of repose, while the other did not, remaining cylindrical. 
When the pile collapsed, its runoff obeyed the scaling 
behavior shown in Fig. [21 

For all aspect ratios below 24, a critical height exists 
above which the pile always collapsed. Additionally, for 
each of these aspect ratios there exists a minimum height 
below which the pile never collapsed. Increasing the par- 
ticle aspect ratio increased these critical heights, and so 
the initial radius does not seem the appropriate quantity 
by which to characterize this behavior. We find, however, 
that these critical heights are the same for all aspect ra- 
tios tested when scaled by the particle length L. This 
is shown in Fig. [3] for acrylic particles. The dashed lines 



FIG. 3: The critical heights that determine whether a pile 
acts as a solid or granular material, scaled by particle length, 
are independent of aspect ratio. The solid and granular re- 
gions are separated by a transition region in which the pile 
has a probability of collapsing, although there are no visible 
differences between piles that collapse and those that do not. 



The probability for a pile to collapse increases as one 
moves through the transition region. We use the upper 
and lower critical heights Hi and Hu to parametrize the 
height, defining an height order parameter that is at 
Hi and 1 at Hu- 



H = 



H-Hi 
H,, — Hi 



(4) 



In Figure [5] we plot the collapse probability as a function 
of this order parameter H for all of our data, encompass- 
ing the complete range of aspect ratios that collapse (4- 
20) and two different particle diameters. Figure |4] shows 
that, when scaled as per Eq. 21 all the data collapse onto 
a straight line between (0,0) (the point at which no piles 
collapse) to (1,1) (where all piles collapse). The data in 
Fig. |4] are grouped by height into bins of scaled width 
0.1, corresponding to the horizontal error bars. The ver- 
tical error bars represent one standard deviation of the 
average collapse probability of all data in that bin. 

Figures[3]and|4]represent the first quantitative descrip- 
tion of the transition region in which the material can 
behave either as a solid or as a granular material. This 
region is unique to large aspect-ratio materials, and ap- 
pears to exist even when the aspect ratio is moderate 
(^ 4). Computational work is underway to determine 
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FIG. 4: The probability of collapse within the transition re- 
gion increases linearly with pile height scaled by the criti- 
cal heights. This behavior is independent of particle length, 
width, and aspect ratio. The normalization by critical heights 
is done so that the origin is the lower transition point and (1,1) 
is the higher transition point. 



what internal characteristics (e.g. distribution of forces 
or percolating force chains) might be responsible for the 
different behavior in superficially similar piles. 

The region of bistability was observed in two differ- 
ent cylinder radii, and the critical heights are normal- 
ized by the cylinder radius. Both cylinders were larger 
than a particle length; it is not clear how piles formed 
in extremely small containers (compared to the parti- 
cle length) would respond. In this case particles could 
not orient horizontally, and this may reduce the stabil- 
ity of any pile formed. Desmond et al. [20| found that 
extremely large containers allowed for a more granular 
response to an intruder, but it is not clear that the same 
behavior exists in this geometry. We saw no decrease in 
pile rigidity in the larger cylinder tested. 



V. CONCLUSIONS 

We have investigated the collapse of granular columns 
consisting of cylindrical rods. Interestingly, when the 
pile collapses, it does so in a manner very similar to that 



of spheres or irregular sand. The runoff distance, when 
scaled by the initial pile radius, initially scales linearly 
with the normalized pile height. This corresponds to the 
partial collapse of the pile, involving only the material at 
the outer edges. The final pile shape is that of a trun- 
cated cone, with outer radius determined by the angle 
of repose of the material. At a critical pile height, also 
determined by the angle of repose, there is a crossover to 
a different scaling, where the runoff grows as the square- 
root of scaled height. This indicates an immobile inner 
cone of material, at the angle of repose, and that the final 
pile height is set by the initial radius. Both the scaling 
behavior, and the location of the cross-over height, are 
quantitatively consistent with earlier work on spheres, 
and so we conchide that, when rodpiles collapse, they 
are quantitatively the same as ordinary sand or spheres. 
Large-aspect ratio particles do display a fundamentally 
different behavior from that of spheres, however, in their 
ability to maintain the shape of their initial container, 
essentially acting as a solid. We have shown that the 
onset of this behavior with increasing aspect ratio is not 
sudden, but manifests in even moderate aspect ratios un- 
der suitable experimental geometries. For a given aspect 
ratio there exist critical pile heights that demarcate the 
solid and granular behaviors, and the transition region 
separating the two. The transition region is character- 
ized by a probability of solid/granular response, with a 
linear transition between the two. The transition heights, 
when scaled by particle length, are independent of parti- 
cle aspect ratio. This behavior persists until the aspect 
ratio is quite large (> 24), at which all piles are essen- 
tially solid. Future research will investigate the internal 
characteristics that determine the pile response, as well 
as the behavior of the critical heights with increasing par- 
ticle aspect ratio and pile radius. 
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